Non-oriented Solutions of the Eikonal Equation 



Mark A. Peletier and Marco Veneroni 
November 24, 2008 

Abstract 

We study a new formulation for the eikonal equation |Vu| = 1 on a bounded sub- 
set of M^. Instead of a vector field Vu, we consider a field P of orthogonal projections 
on 1-dimensional subspaces, with divP € L^. Wc prove existence and uniqueness for 
solutions of the equation Pdiv-P = 0. We give a geometric description, comparable 
with the classical case, and we prove that such solutions exist only if the domain is 
a tubular neighbourhood of a regular closed curve. The idea of the proof is to apply 
a generalized method of characteristics introduced in [8] to a suitable vector field m 
satisfying P = m. 

This formulation provides a useful approach to the analysis of stripe patterns. 
It is specifically suited to systems where the physical properties of the pattern are 
invariant under rotation over 180 degrees, such as systems of block copolymers or 
liquid crystals. 

AMS CI. 35L65, 35B65. 

1 Introduction 

1.1 Stripe patterns and the eikonal equation 

Many pattern-forming systems produce parallel stripes, sometimes straight, sometimes 
curved. In geology, for instance, 'parallel folding' refers to the folding of layers of rock 
in a manner that preserves the layer thickness but allows for curving of the layers [3]. 
In a different context, the convection rolls of the Rayleigh-Benard experiment produce 
striped patterns that may also be either straight or curved (see e.g. [2]). The system 
that suggested the work of this paper is a third example: in [9] we investigated striped 
patterns that arise in the modelling of block copolymer melts. 

Block copolymers consist of two covalently bonded, mutually repelling parts ('blocks'). 
At sufficiently low temperature the repelhng forces lead to patterns with a length scale 
that is related to the length of single polymers. We recently studied the behaviour of an 
energy that describes such systems, and investigated a limit process in which the stripe 
width tends to zero [9]. In that limit the stripes not only become thin, but also uniform in 



width, and the stripe pattern comes to resemble the level sets of a solution of the eikonal 
equation. The rigorous version of this statement, in the form of a Gamma-convergence 

result, gives rise to a new formulation of the eikonal equation, in which the directionality 
of the stripes is represented by line fields rather than by vector fields. Before stating this 
formulation mathematically we first describe it in heuristic terms. 

The eikonal equation has its origin in models of wave propagation, where the equation 
describes the position of a wave front at different times t. For a homogeneous and 
isotropic medium, in which the wave velocity is constant, the equation can be written in 
the form 

\Vu\ = 1. (1) 

The wave front at time t is given by the level set {x : u{x) = t}, and the function u has 
the interpretation of the time needed for a wave to arrive at the point x. 

A feature of the eikonal equation is that the fronts at different times are parallel, 
provided no singularities occur. In this sense the equation is a natural candidate for the 
description of other processes that involve parallellism, such as the stripe-forming systems 
mentioned above. However, a major difference between the stripe-forming systems and 
the wave-front model is that the wave front has a natural directionality associated with it: 
of the two directions normal to a front, one is 'forward in time' and the other 'backward'. 
This distinction also is visible in the notion of viscosity solution for (1) (see e.g. [4]). 

The stripe patterns, on the other hand, have no inherent distinction between the two 
normal directions. As a consequence a vector representation of a stripe pattern may have 
singularities that have no physical counterpart. Figure 1 (left) shows an example of this. 




Figure 1: Stripe patterns can be represented by vectors (left) or by unoriented line fields 
(right). Both representations have a vortex singularity at the two ends; but the vector 
representation also contains a jump singularity along the connecting line. Note that the 
regularity restrictions of this paper exclude both types of singularity, however. 

In our work on block copolymers that was mentioned above [9], the limiting stripe 
problem is formulated in terms of unsigned vector fields, or line fields, which capture 
direction only up to a sign (Figure 1 (right)). These stripe patterns are required to 
satisfy what is essentially the eikonal equation for line fields (see (2) below). In this 
paper we investigate this new formulation of the eikonal equation, study the properties 
of its solutions, and compare it to various existing formulations. We only consider the 
case of two space dimensions. 
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1.2 The eikonal equation in terms of line fields 

We now turn to the work of this paper. A natural mathematical object for the represen- 
tation of line fields is a projection. For the purposes of this paper, we define a projection 
to be a matrix P that can be written in terms of a unit vector m as P = m ® m. 
Such a projection matrix has a range and a kernel that are both one-dimensional, and if 
necessary one can identify a projection P with its range, i.e. with the one-dimensional 
subspace of M? onto which it projects. Note that the independence of the sign of m, 
that is the unsigned nature of a projection, can be directly recognized in the formula 
P = m® m. 

Wc define div P as the vector- valued function whose i-th component is given by 
(divP)i := Yl^j=i'^xjPij- We consider the following problem. Let Q. be an open subset 
of M^. Find P G L°°(J^; M2x2^ g^^j^ ^-^^^ 

P"^ = P a.e. in Q, (2a) 

rank(P) = 1 a.e. in 9., (2b) 

P is symmetric a.e. in $7, (2c) 

div P G L^(m2; M^) (extended to outside 9.), (2d) 

P divP = a.e. in 9. (2e) 

The first three equations encode the property that P{x) is a projection, in the sense 
above, at almost every x. Equation (2d) is both a regularity requirement and a boundary 
condition, and the choice of the exponent 2 has a precise explanation, as we show below. 

Given the regularity provided by (2d), the final condition (2e) is the eikonal equation 
itself, as a calculation for a smooth unit-length vector field m{x) shows: 

= P div P = m{m ■ {m div m + Vm • m)) = m div m -|- m{m ■ Vm ■ m) = m div m, (3) 

where the final equality follows from differentiating the identity |mp = 1. A solution 
vector field m therefore is divergence- free, implying that its rotation over 90 degrees 
is a gradient Vu; from |m| = 1 it follows that IV^I = 1 (see Section 5.2). This little 
calculation also shows that the interpretation of m in P = m (g) m is that of the stripe 
direction; P projects along the normal onto the tangent to a stripe. 

The sense of property (2d) is that the divergence of P (extended to outside 0,), in 
the sense of distributions in M^, is an L^(M^) function, i.e. there exists C > such that 
for any test function ip e C° 



< C\Ml2(^2). (4) 



/ P{x) : Vcp{x) dx 
Since the trace Pn satisfies the equality 

- [ P ■.V(pdx= [ dWP-(pdx- [ (Pn) ■ (pdS, (5) 
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condition (2d) implies 

Pn = in the sense of traces on d^}. 
In order to see this, note that by (4) and (5) there exists C > such that 



(6) 



div P ■ ipdx 



an 



{Pn)-ipdS 



< C\\(p\\l2 



(R2)' 



Choosing a sequence {tpk} such that ipk ^ in L^(M^) and (pk — Pn on L'^{dCl), we 
conclude (6). 

The exponent 2 in (2d) is critical in the following sense. Obvious possibilities for 
singularities in a line field are jump discontinuities ('grain boundaries') and target pat- 
terns (see Figure 2). At a grain boundary the jump in P causes divP to have a line 





(a) grain 
boundary 



(b) target and U-turn patterns 




directional 
variation 



Figure 2: Canonical types of stripe variation in two dimensions. Types (a) and (b) are 
excluded by (2d). 

singularity, comparable to the one-dimensional Hausdorff measure; condition (2d) clearly 
excludes that possibility. For a target pattern the curvature k of the stripes scales as 1/r, 
where r is the distance to the center; then f is locally finite for p < 2, and diverges 
logarithmically for p = 2. The cases p < 2 and p > 2 therefore distinguish between 
whether target patterns are admissible {p < 2) or not. 



1.3 Results 

The first main contribution of this paper is to show that problem (2) indeed resembles 
the eikonal equation. In the heuristic discussion above, we argued that the field P has 
to be locally parallel, to be parallel to the boundary of fi, and to avoid line and vortex 
singularities. The two first statements are formalized in the following theorem. 

Theorem 1. Let be an open, bounded, and connected subset o/M^ with boundary, 
and let P be a solution of (2). 

1. Let xq be a Lebesgue point of P in O, let x € 11, and let L be the line segment 
connecting xq with x. Assume that L C fi. // P{xo) • {x — xq) = 0, then P{y) = 
P{xo) for Jif^ -almost every y E L. 
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2. P ■ n = a.e. on dft. 

These two statements are meaningful since 

We recall that x is called a Lebesgue point of m if 

lim / \m(x) — m(y)\dy = 0. 

A standard result yields that a.e. x G O is a Lebesgue point for m (see e.g. [7, Section 
1.7]). 

The second main result is to show that the restrictions on P are so rigid that the 
mere existence of a solution provides a strong characterization of the geometry of the 
domain 17: 

Theorem 2. Let O he an open, bounded, and connected subset of M? with boundary. 
Then there exists a solution of (2) if and only if is a tubular domain. In that case the 
solution is unique. 

A tubular domain is a domain in that can be written as 

o = r + 5(o,(5), 

where F is a closed curve in with continuous and bounded curvature k, < 5 < ||av||^^, 
and -6(0, 5) is the open ball of center and radius 5. 




Figure 3: If tangent directions propagate normal to themselves (Theorem 1.1), and if in 
addition the boundary is a tangent direction (6), then the domain is tubular (Theorem 2). 

The reason why Theorem 2 is true can heuristically be recognized in a simple picture. 
Figure 3 shows two sections of dVt with a normal line that connects them. By the first 
assertion of Theorem 1, the stripe tangents are orthogonal to this normal line; by the 
second, this normal line is orthogonal to the two boundary segments, implying that the 
two segments have the same tangent. Therefore the length of the connecting normal line 
is constant, and as it moves it sweeps out a full tubular neighbourhood. 
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1.4 Discussion 



The work of this paper represents a first step in the analysis of this projection- valued 
eikonal equation. While the main results are still lacking in various ways — ^which we 
discuss in more detail below — the main point of this paper is to show that this projection- 
valued formulation is a useful alternative to the usual vector-based formulation. 

To start with, our Theorems 1 and 2 show that solutions of (2) behave much like we 
expect from the eikonal equation, in the sense that directional information is preserved 
in the normal direction. Theorem 2 makes this property even more explicit, by showing 
that a full tube, or bunch, of parallel 'stripes' can be identified. 

However, it is the differences with the vector-valued eikonal equation that are the 
most interesting. Figure 1 shows how this formulation can be a better representation of 
the physical reality than the vector-based form. On the left, the vector field has a jump 
discontinuity along the center line, while on the right the projection is continuous along 
that line. Depending on the underlying model, this singularity may have a physical 
counterpart, or may be a spurious consequence of the vector-based description. For 
the wave-propagation model the singularity is very real; for striped-pattern systems it 
typically is not. A projection- valued formulation therefore provides an alternative to the 
Riemann-surface approach that is sometimes used [6]. For this distinction to have any 
consequence, however, solutions with less regularity than the divP G of this paper 
are to be considered. 

We now comment in more detail on our method of proof. The proof of the properties 
that we give in this paper relies on a reduction of the projection-valued formulation 
to a vector-based formulation. This reduction is achieved by the Ball-Zarnescu lemma 
(Lemma 3) , which requires div P E L^; for less regularity the existence of a lifting may 
not hold, as the example of the U-turn pattern (Figure 2b) shows. 

The dependence of the proof on a vector-based representation is awkward in various 
ways. To start with, the condition divP € required for the lifting is much stronger 
than the conditions (14-16) that Jabin, Otto, and Perthame require for their results [8]. 
It also has the effect of excluding all singularities, as we already remarked. It would be 
interesting to prove properties such as those of Theorems 1 and 2 by methods that do 
not rely on this lifting. 

We would hope that such an intrinsic projection-based proof could also be generalized 
to the study of target patterns and U-turns, and eventually of grain boundaries. These 
will require increasingly weak regularity requirements: target patterns may exist for 
divP G LP with p <2, and for a line discontinuity, such as a grain boundary, divP will 
be a measure. 
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1.5 Overview of the paper 



The structure of the rest of the paper is as follows. We start, in Section 2, by showing 
that P G (part 3 of Theorem 1). In Section 3 we construct a lifting, i.e. a vector field 
m such that P = m®m\ there we also show part 2 of Theorem 1. In Section 4 we show 
that the lifting m satisfies the conditions of [8] , and prove the final part of Theorem 1 . 
In Section 5 we combine all the earlier results to prove Theorem 2. 



2 Regularity of P 

The regularity statement 3 of Theorem 1 follows from simple manipulation. By (2a-2c) 
P is a symmetric matrix which we can write as 

a b 
b c 



and whose elements satisfy 



a,b,ce L°°{n), 0<a,b,c<l, (7) 
b^ = ac, a + c=l. (8) 



(9) 



Denote by Oj := S^r-a (idem for bi,Ci) for i = 1,2. By (2d) there exist /, 5 G L^(n) such 
that 

ai + ^2 = / 
bi + C2 = 9. 

From (8) we deduce by differentiation 

bbi= f ^ - « j ai> 



and by using this and (7-9) in (2e) we have 



= a{ai + b2) + b{bi + C2) 

= a{ai + 62) + ~ '^^ + ^"^2 

= ^0,1 + 062 — ba2 

= ^ai + a{f - ai) - ba2. 
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Similarly 



= 6(ai + 62) + c(6i + C2) 

= bai+ (I- — c] C2 + c{bi + C2) 



= bai - ^02 + (1 - a)bi 

= 6ai - ^02 + (1 - a)(cj' + 02). 

Therefore ai , 02 solve the Hnear system 

^-a -b \ f 0-1 \ ^ ( -«/ 
b l-a J \ a2 J \ {a-l)g 

Since the determinant of the coefRcient matrix is — a)^ + b^ = ^ > 0, the system is 
nondegenerate, and it follows from a,b,c & L°°(Q) and f,g ^ that ai,a2 G L^(ri). 

By the relations (8) and (9) we find ci, C2 G L^(0) and 61, 62 G L'^{^). Therefore 

PijeH\n), i,j = 1,2. 

□ 



3 Lifting to a vector representation 

The remaining parts of both theorems are proved by using a vectorial representation 
of P. We say that P is orientable if there exists a vector field m G H^{Cl; S^) such that 

P = m® m a.e. on fi. (10) 

In this case m is called a lifting of P. Similarly, the trace of P on the boundary is defined 
to be orientable if there exists a vector field m G H^/^{dfl; S^), such that 

Tr(P)|gf2 = m m a.e. on 5$!. (11) 

The following lemma by Ball and Zarnescu establishes a link between orientability in 
the bulk and on the boundary. We formulate it in the language of this paper. 

Lemma 3 ([1]). Let J7 be an open, bounded, connected subset o/M^, with boundary. 
If P satisfies (2) and has the additional regularity P G i?^(r2; M^^^), then P is orientable 
if and only if its trace on dO, is orientable. 

Note that the boundary dO, has only a finite number > 1 of connected components. 
Indeed, let A be a set of indexes, and let 7a, for A G A, be a connected component of 
dCl, which disconnects into a bounded set Bx and an unbounded set U\. Since CI is 
connected. 
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• there is a unique index A such that B^Ci^ ^ $ and 

• Bx, n = if Ai ^ A2 and Ai, A2 / A. 

Since dQ is C^, its curvature is bounded by some constant 6 and therefore every Bx must 

contain at least one ball of radius 1/S. Since dfl is bounded, UagAAt^a-^a is bounded 
and only a finite number of such balls can fit into it without overlapping. We conclude 
that #A < 00. 

We now construct a lifting m of P on the boundary. Let dCl = dCf U • • • U be a 

decomposition of the boundary into connected components, and let aj : [0,Lj] dQ.\ 
be arclength parameterizations of dO^K Then, owing to (6), the vectors defined 
by fn^{aj{s)) := Oi'j{s), j = . . . N, satisfy (11). Therefore, by Lemma 3 we obtain that 
P is orient able. 

Thus we have proved 

Lemma 4. There exists a vector field m € H^{fl,'R'^) such that P = m® m. 

This proves part 3 of Theorem 1. Part 2 then follows from the calculation (5). Note 
that Pn = on the boundary implies 

m - n = a.e. on dCl. (12) 



4 Properties of m 

We use some concepts and notation from [8]. We define the functions x ^-iid X by 

/ / / / 1 if m{x) • C > , , 

X{x,0--=x{m{x),0:=[^ if m(a;) • ^ < 0. ^'^^ 

Lemma 5. The vector field m satisfies 

\m{x)\ = 1 for a.e. xinCl, (14) 
div m = distrihutionally in M^, (15) 
^•Vx(',C) = distrihutionally in for all in . (16) 

Proof. Property (14) follows from remarking that 

\m\^ = m- P- m = m-P^-m = \mf' . 

The computation (3) yields (15). The proof of (16) is an application of some ideas 
introduced in [5] and further developed in [8]. We use the notion of entropy introduced 
in [5, Definition 2.1]. A vector field $ G Cq°(M^,M^) is called an entropy if 

z-V^z^ = Q forallzGM^ and $(0) = 0, V$(0) = 0. (17) 
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In [5] and [8] it was shown that for every entropy $ there exists a vector field ^ G 
C^(M^M2) g^^^ ^ function a G C^(M^M) such that for any m e H^{n,M.'^) 

div($(m)) = ^'(m) • V(l — |m|^) + a{m) divm a.e. in 

As proved in [5, Lemma 5], for any fixed ^ G 5^, the function 

R''3z^x{z,0C (18) 

(where x is defined in (13)) is the pointwise limit of a sequence {$„}neN of entropies in 
the sense of (17). Properties (14-15) and the approximation (18) then yield equation (16). 

□ 

The properties of m stated in Lemma 5 allow us to apply Proposition 3.2 in [8], which 
reads 

Lemma 6 ([8]). Let m satisfy (I4-I6). Let G O 6e a Lebesgue point for m in Q, and 
let L G O, be a straight line segment containing xq. Then 

m{x) ■ m'^{xo) = for J^^-a.e. x E L. 

In terms of P = m(8)m this statement reduces to part 1 of Theorem 1. This concludes 
the proof of Theorem 1. 

5 Proof of Theorem 2 

The statement of equivalence in Theorem 2 contains a trivial and a non-trivial part. The 

non-trivial part is to show that existence of a solution P implies that 0, is tubular. The 
trivial part is to construct a solution P, if one assumes that 0, is tubular. We prove the 
non-trivial part first, since the calculations will be useful in the second part. Part 1 of 
Theorem 1 will appear as an intermediate result in the proof of Theorem 2. 

5.1 A first chciracterization of Q 

Let A be the class of sets CI such that 

O C is open, bounded, connected, C^, and 3y,z e dCl such that the nor- 
mal hnes issued from y, z are different and intersect in CI before crossing dCl. 

Then Theorem 1.2 in [8] states that ii CI e A and m satisfies (14), (15), and (16), 
then n is a disk and m is a vortex, i.e. there exist Xq e CI and a = ±1 such that 

(x — Xq) 

mix) = a— j- for a.e. x eCI. (19) 

\x — Xo\ 

(Here and below, is the counterclockwise rotation by 90 degrees of the vector v.) A 
direct computation shows that if P = m (8> m with m given by (19) then divP ^ L'^{Cl). 
It follows that CI ^ A, and therefore either the normal lines issued from any two points 
on the boundary are identical or they have no intersection in CI. 
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5.2 Construction of a potential 

Let 7 : [0, L] — >■ dil^ be a arclcngth parameterization of d^l^ and let n be the exterior 
unit normal vector on dCl. By property (15) the vector field m"*" satisfies 

curl m"*" = 0, distributionally in M^. 

Since the vector field m-*- is orthogonal to the boundary (see (12)), by Green's theorem 
the integral of m-^ along any one-dimensional closed curve in CI is zero. Therefore there 
exists a potential (p G H^i^l) which satisfies 

V(j){x) = m^{x). (20) 

We can assume that cf) is continuous on CI. Prom (12) we deduce 

• n"*" = — m • n = a.e. on dCl, (21) 

therefore (f) is constant on every connected component of dCl: i.e. there exist Cj G M, 
j = 0, . . . , N such that 

(/)(x) = Cj \/x e dCl^, j = l,...,N. (22) 

It is not restrictive to assume cj > 0, cq = 0. Moreover for x G dCl'^ and T > we define 
the segments 

Lx,T ■= {x-tnix) ■.te[0,T]}, 

which intersect dCl orthogonally by (21) (wherever they intersect), and the function 
T{x) : dn° M, 

T{x) := max{T > : L^^t C U}. 
Note that it follows from this definition that 

X - T{x)n{x) e dn. (23) 

5.3 O is tubuletr 

By [7, Theorem 2, section 5.3], there exists a set A/" C dCl such that J^^{M) = and 
every x G dQ\M is a Lebesgue point for m with respect to the two-dimensional Lebesgue 
measure. 

We now proceed in five steps. 

«^ As a direct consequence of Lemma 6, 

V.T G dCl,'^\M, for Jf^-a.e. y G L^^j^(^^^, m{y) = ±m(x); 
a) We claim, in addition, that 

Va; G dCl^\J\f, for J^-'^-a.e. y G ija;,T(a;)) "miv) = m{x); 
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otherwise, following the arguments in the proof of [8, Theorem 1.1], it can be seen that 
a point where ?n.(y) jumps from m(x) to —m{x) would be a vortex point for m; as we 
remarked above, that would imply that divP L^, a contradiction. 

in) By (20) and step ii), for each xq £ 9^2*^ \ AA, (f) is linear on L^^^^^^^-^ with slope 1. 
iv) By (23) and step in), for J^^-a.e. x G dQ^ there exists a j S {1, . . . , A^} such 
that 

c j — = (/)(x — r(x) n(x)) — 4){x) = T(x). 
Therefore there exists a step function T : — > {ci, . . . , cat} such that 

T = f, jr^-a.e. on 517°. 

Let X be a discontinuity point for T. 




Figure 4: For example, T in x jumps between the values ci and C2- 

We claim that there exists x € i^x,T(s) ^ such that n{x) ■ n(x) = 0; otherwise, the 
regularity of dQ and the continuity of n would make T continuous in x. This contradicts 
^ ^ A, and we conclude that T is continuous and T{x) = c. Finally, the image of d^l^ 
under the continuous function x i— > x + T(x) n{x) is connected, and therefore dO, consists 
of two connected components at constant distance c: dQ = (?J7°UOr2^ and O = r+S(0, 6) 
for 6 = c/2 and F := {x+ |n(x), x £ 917°}. This concludes the proof of one half of Theo- 
rem 2. □ 

5.4 Proof of the converse 

Let O be a tubular domain of width 6 and let 50 = 50° U 90^ be the decomposition 
of the boundary into connected components. Let 7 : [0, L] — > 90° be a arclength 
parameterization of 50°, then 

^ : [0,L] X (0,25) ^ O 
^'(s,t) := 7(5) - tn{s) 

is a parameterization of O. Then the line field P defined as 

P{^{s,t)):=j'{s)CSj'{s) (24) 



12 



is a solution and satisfies P G C^(0,M^^^). Since by step ii) above m is uniquely deter- 
mined a.e. on CI by its trace on dCl, i.e. by the tangents to dCl, we conclude that P is 
uniquely determined a.e. on CI and (24) is the unique solution. 
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